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Abstract 

A Quintom universe with an equation-of-state crossing the cosmological constant boundary can 
provide a bouncing solution dubbed the Quintom Bounce and thus resolve the Big Bang singularity. 
In this paper, we investigate the cosmological perturbations of the Quintom Bounce both analytically 
and numerically. We find that the fluctuations in the dominant mode in the post-bounce expanding 
phase couple to the growing mode of the perturbations in the pre-bounce contracting phase. 
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1 Introduction 

For years, it has been suggested that bouncing universe scenarios could provide a solution to the problem 
of the initial singularity in the evolution of the cosmos, a problem which afflicts both Standard Big 
Bang cosmology and inflationary universe models [U [51 [31 H] 0. In the literature there have been a 
lot of efforts towards constructing bouncing universes, for instance the Pre-Big-Bang scenario [7] or the 
cyclic/Ekpyrotic scenario [5]. Also, in loop quantum cosmology the evolution might be singularity-free [S]. 
In Refs. [ini [HI [12] and [131 HI] models were considered in which the gravitational action is modified by 
adding high order terms. Since, in general, in these models the equations for cosmological perturbations 
are very involved, it is difficult to make definitive conclusions about the evolution of fiuctuations, and the 
knowledge of the evolution is crucial in order to compare the predictions of the model with observations. 

In some previous work [15j . a subset of the present authors considered a bouncing universe model 
("Quintom Bounce") obtained within the standard 4-dimensional Friedmann-Robertson- Walker (FRW) 
framework, but making use of Quintom matter [16] . In this model, it can be shown in [15j that the 
null energy condition (NEC) is violated for a short period around the bounce point. Moreover, after the 
bounce the equation-of-state (EoS) of the matter content w in our model is able to transit from w < —I 
to w > — 1 and then connect with normal expanding history. 

A Quintom model jl6| was initially proposed to obtain a model of dark energy with an EoS parameter 
w which satisfies w > —1 in the past and w < —1 at present. This model is mildly favored by the 
current observational data fitting [17| . The Quintom model is a dynamical scenario of dark energy with 
a salient feature that its EoS can smoothly cross over the cosmological constant barrier w = —1. It 
is not easy to construct a consistent Quintom model. A no-go theorem proven in Ref. (TH] (also see 
Refs. [ini [ini [201 mi im us HU) forbids a traditional scalar field model with Lagrangian of the general 
form £ = C{4>, d^cjjd^ff)) from having its EoS cross the cosmological constant boundary. Therefore, it is 
necessary to add extra degrees of freedom with un-conventional features to the conventional single field 
theory if we expect to realize a viable Quintom model in the framework of Einstein's gravity theory. The 
simplest Quintom model involves two scalars with one being Quintessence-like and another Phantom-like 
[TSl [IS] . This model has been studied in detail in . In recent years there have been a lot of theoretical 
studies of Quintom models. For example, motivated from string theory, the authors of Ref. [2^ realized 
a Quintom scenario by considering the non-perturbative effects of a DBI model. Moreover, there are 
models which involve a single scalar with higher derivative terms in the action 28J, models with vector 
field [21], making use of an extended theory of gravity [30], of ideas from non-local string field theory 
theory [3T], and others (see e.g. [311 [331 [33] ) . Because of its brand new properties, Quintom models give 
many unexpected predictions. Most important for the subject of the present paper is that a universe 
dominated by Quintom matter can provide a bouncing cosmology which allows us to avoid the problem 
of the initial singularity. This scenario, developed in [15j, is called the Quintom Bounce. In this scenario, 
the EoS crosses the cosmological constant boundary twice around the bounce point. 

In order to be able to compare the predictions of a Quintom bounce with observations, it is crucial 
to investigate the evolution of cosmological perturbations. In singular bouncing models the evolution 
of fiuctuations is not under control classically. The fluctuations diverge at the bounce point, firstly 
invalidating the applicability of the theory of linear cosmological perturbations, and secondly making it 

II An initial singularity occurs in inflationary models if scalar fields are used to generate the accelerated expansion of 
space [5]. If inflation arises as a consequence of modifications of the gravitational part of the action, as in Starobinsky's 
original proposal [6], it is possible that the evolution is singularity- free. 
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impossible to reliably connect the fluctuations in the contracting phase with those in the expanding phase. 
In a Quintom bounce model the evolution of fluctuations is well-behaved. Our first main result is that 
cosmological perturbations can evolve smoothly from the contracting to the expanding phase. We show 
that the commonly used variable the metric perturbation in longitudinal gauge, passes smoothly both 
through the bounce point and through the points where the cosmological constant boundary of the EoS is 
crossed. Interestingly, we find that our result for the transfer of fluctuations through the bounce possesses 
features found both in some previous analyses of fluctuations in non-singular cosmologies, and features 
of the evolution in singular bouncing models. For small comoving wave numbers one finds an evolution 
similar to what occurs in singular bouncing cosmologies (the growing mode of $ in the contracting phase 
couples only to the decaying mode in the expanding phase), whereas the transfer of the perturbations is 
very different in the region of large comoving wave numbers (the growing mode in the contracting phase 
couples to the dominant mode in the phase of expansion). The physical reason is easy to see: small 
k modes enter the sub-Hubble region very close to the bounce point and have no time to complete a 
single oscillation, whereas large k modes undergo many oscillations which are in the bounce phase at sub- 
Hubble scales. Thus, whereas the large k modes feel the resolution of the singularity, the small k modes 
evolve as if there still were an abrupt transition between the contracting to the expanding phase. We will 
follow the fiuctuations both numerically and analytically. In the analytical analysis, we apply matching 
conditions away from the bounce point, in regions where both the background and the perturbations can 
be matched consistently via the Deruelle-Mukhanov [32] (initially derived in ) conditions Ell- 

This paper is organized as follows. In Section 2, we review the equations of motion of double-field 
Quintom and derive the perturbation equations for this model. In Section 3, we study the behavior of a 
specific Quintom Bounce model systematically. In Subsection 3.1, we determine the complete evolution 
of the background using numerical integration. In Subsection 3.2, we investigate the perturbations, first 
analytically by solving the equations in separate periods by means of approximations valid during those 
periods, and matching the solutions at the transition points, and finally numerically. The numerical 
results support our analytical calculations consistently for all times. Section 4 contains discussion and 
conclusions. 



2 Review of the Quintom Bounce in a Double-Field Model 

2.1 Equations of Motion of Double-Field Quintom 

To start, we take a Quintom model consisting of two fields with the Lagrangian 

where the signature of the metric is (+, — , — , — ). Here the field (j) has a canonical kinetic term, but -0 
has a kinetic term with the opposite sign and thus plays a role of a ghost field. In the framework of a 
flat Friedmann-Robertson- Walker (FRW) universe, the metric is given by ds^ — dt^ — a? {t)dx'^ dx^ . By 
varying the corresponding matter action, we easily obtain the following expressions for the energy density 
p and pressure p of this model, 

p^\4>^- + ^('^) + ^(V-) , 'p = \^^- ^V-' - v{<i>) - w{^) , (2) 

**For a discussion of the dangers in applying the matching conditions at the transition point between contraction and 
expansion see |37| . 
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and the background equations of motion arc given by 



H^ = ^[U^-^4'^ + V + W], (3) 



3 1-2^ 2 

j,2 J,2\ 



H = -47rG(0^ - , (4) 
<^ + 3770 + V;^ = , (5) 
i/i + 3Hi; -W.^ = , (6) 

where a{t) is the scale factor, H = ^ is the Hubble parameter and the dot denotes the derivative with 
respect to the cosmic time t. 

It is well known that this model can realize the EoS it; of a universe which crosses the cosmological 
constant boundary. Namely, the expression for w is given by 

w = —I + —. —. . (7) 

1 02 _ 1^2 + + 

If (p^ — tp'^ is positive, w is larger than —1 and the universe lies in a Quintessence- like phase; if 0^ — -0^ 
is negative, w is less than —1 and correspondingly the universe is Phantom- like. In recent work, it was 
found that the Quintom model can provide a bouncing universe scenario (15). If we work in the context of 
a spatially flat four-dimensional background FRW metric, there must be a period when the Null Energy 
Condition (NEC) is violated if we are to obtain a smooth transition from a contracting universe into an 
expanding phase. In this case, we need a kind of matter which admits an EoS parameter which is less 
than —1, but only around the bounce. Neither regular not Phantom matter alone can achieve a transition 
in the EoS parameter through the cosmological constant boundary. Therefore, a Quintom model is the 
only possible solution to resolve this difficulty. 

2.2 Equations of Motion of Perturbations 

Now let us consider linear perturbations of the metric. The longitudinal (conformal-Newtonian) gauge 
metric perturbations are given by 

ds^ = a^{ri)[{l + 2^)d'n^ - (1 - 2^)dx'dx'] , (8) 

where we introduced the comoving time rj defined by drj — dt/a for convenience. The fields $ and ^' 
depend on space and time and contain the information about the scalar metric fluctuations, the degrees 
of freedom in the metric which couple to matter fluctuations (see e.g. [38) for a comprehensive survey 
of the theory of cosmological perturbations and ^\ for a recent overview). We will not discuss vector 
fluctuations since they are not induced by scalar field matter. To linear order, tensor metric fluctuations 
(gravitational waves) do not couple to matter and hence we will not discuss them here. By expanding 
the Einstein equations to first order, we obtain the following perturbation equations 

V^* - 37i:(^'' -t- n<i>) = AnGa^Sp , (9) 
+ = -AirGaSq , (10) 

+ 37i:$' + (27i' + = AnGa^Sp , (11) 

where Ti. = da/adrj is the comoving Hubble parameter and the prime denotes the derivative with respect 
to the comoving time. Here dp, Sq and dp represent the perturbations of density, momentum and pressure, 
respectively. The first equation is the perturbed '00' equation, the second the perturbed 'Oi' equation, and 
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the third is the diagonal 'ii' equation. For scalar field matter there is no anisotropic stress to linear order 
in the matter fluctuations, and hence it follows from the ofii'-diagonal 'ij' perturbed Einstein equations 
that $ = 

From tlic! conservation of the stress tensor, we also obtain the analog of the continuity equation for 
the matter fluctuations: 

V + m{5p + 5p) - i V^^g = 3(p + . (12) 

In the double-field Quintom model, the scalar fields are perturbed: (j) ^ (p + 5(p and ip ^ ip + 6ip. 
To linear order in the scalar field fluctuations and dtp, the perturbations of energy, momentum and 
pressure can be expressed by 

5p = ^(l>'{5(p' - 0'$) + V^S(l> - ^V' W - Z^) + W,^6ij , (13) 

Sq = ^[-^'S(l) + ^'S'tp] , (14) 

6p = -^,p'{6(P' -<j3'^)-V4,6(b-^^^;'{5^' -^P'^)-W,^6^p , (15) 

respectively. The linear terms in the equations obtained by varying the matter action with respect to the 
two matter flelds, we obtain the equations of motion for these scalar field perturbations, 

S(f)" = -2nS(j)' - k'^Scj) - a'^V^^Scp - 2a'^V^^ + 4^'$' , (16) 
6ip" = -2H6ip' - k'^Sip + a^W^^^Sip + 2a'^W,^^ + 4^^'$' . (17) 



3 A Concrete Example of a Quintom Bounce 

In order to make the analysis more quantitative, in this section we would like to take a specific example 
of a Quintom model. Moreover, our Quintom model should not only provide a bouncing solution of 
the universe, but also realize a scenario where the EoS crosses the cosmological constant boundary and 
yields a normal (i.e. non-Phantom-like) evolution of the universe at very early and very late times. 
Consequently, we demand that the ghost field ip dominates only close to the bounce point. It must be 
able to decay after the bounce. 

3 . 1 Background 

We will assume vanishing potential for the ghost field, whose background energy density thus consists 
of kinetic energy exclusively. From the equation of motion of the ghost field it follows that ip evolves 
proportional to a~^. Thus, the kinetic energy density of the ghost field is proportional to a~^. Therefore, 
the contribution of the ghost field tjj will dominate at the bounce point, but will decay quickly before and 
afterwards. This is just what we need in order to achieve both a bounce and connect to the observed 
universe. 

We thus take the Lagrangian of Quintom model to be 

jC = ^d^m - \d^i>d^i^ - V{<f>) , (18) 
and, to be speciflc, we choose a simple form of the potential 
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Figure 1: Plot of the evolution of the Hubble parameter H in the model p^ . In the numerical calculation 
we choose the initial values of parameters as: cj) — —5.6 x 10"^, = 2.56 x 10^, V — 4.62 x 10^^^, m — 
1.414 X 10^^. The initial time was chosen to be t = —500. Note that in this and the following figures 
with results from numerical simulations, all masses are expressed in units of 10~^Mp;. 




Figure 2: Plot of the evolution of the EoS parameter w in the model (fT5|) . The initial values of parameters 
are the same as in Fig. [TJ 
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In this model, the evolution of the universe is symmetric about the bounce point. In the region long 
before the bounce point, we choose initial conditions for which the energy density in the field cj) dominates 
by a large factor. Thus, in the initial phase of evolution (Period 1) the field (p ^nd hence also the EoS 
of the universe will oscillate about w = and the average state is similar to that in a matter dominated 
period. Since the universe is contracting, it is heating up. Once the amplitude of the field oscillation 
gets sufficiently large and the absolute value of the Hubble expansion parameter increases sufficiently, 
the universe will enter a nearly de-Sitter contracting phase (Period 2) with w ~ — 1. This transition is 
the time reversal of the end of the slow-rolling period in scalar field inflation. During this second period, 
the period of quasi-de-Sitter contraction, the field (/) climbs slowly up its potential. During both Periods 
1 and 2, but in particular in Period 2, the contribution of -0 to the energy density is increasing. We take 
initial conditions such that at the transition point between Periods 1 and 2 the contribution of ip is still 
subdominant. Due to the exponential decrease of the scale factor during Period 2, the contribution of -0 
rises exponentially and soon starts to dominate. At this point, the EoS parameter w crosses the barrier 
w = —1. This triggers the onset of Period 3, the bouncing phase. The parameter w approaches negative 
infinity at which time the bounce occurs. After the bounce, the time reverse of the evolution before the 
bounce occurs. The universe again enters a quasi-de-Sitter phase (Period 4), but this time the universe 
is expanding nearly exponentially, i.e., a period of inflation takes place. During this period, (j) is rolling 
slowly but the ghost field decays rapidly, and the EoS returns to a non-Phantom- like value w > —1. Once 
(j) decreases below a critical value (the same critical value which signifies the end of inflation in an m^0^ 
inflation model), it will start to oscillate about the minimum of its potential, and the EoS starts to once 
again oscillate around w = 0. 

In Figs. [T]and[2]we plot the evolution of the Hubble parameter and of the EoS, calculated numerically. 
In this calculation, we normalized the dimensional parameters such as m, _ff , <f> and ■0 by dividing by a 
mass scale M which we took to be 10~^Mp/ where the Planck mass is determined in terms of Newton's 
gravitational constant by Mpi = 1/^/G. From Fig. [21 we can see that in this model the universe undergoes 
a heating phase, a slow-climb-contracting (SCC) phase, a bounce, a slow-roll-expanding (SRE) phase, 
and finally a phase of cooling. The transitions between the phases are smooth. 

Note that the background evolution obtained above is rather generic and does not depend sensitively 
on the details of the initial conditions. The reason is the following: In our model we have introduced 
a Phantom field which makes it rather trivial to obtain a bouncing solution. Since we have chosen the 
potential of the Quintom field to vanish, i.e. W(ip) — 0, the Phantom field evolves as oc which is 
similar to matter with its EoS parameter equal to 1. Thus, we can conclude that a bouncing solution 
is obtained if the absolute value of the Phantom energy density in the contracting universe grows faster 
than that of the normal matter. This will occur if the EoS of the normal matter is smaller than l.In 
our concrete model we have chosen normal matter to be a canonical field (p with a quadratic potential 
V{(j)) = m^0^/2, which leads to an EoS parameter smaller than 1. If the initial conditions arc imposed 
during the period when the Phantom part is negligible, the parameters of the canonical fleld such as the 
mass m and its initial conditions only affect how long the universe is staying in the heating phase when 
w is oscillating around and in the SCC phase when (w ~ —1)- 
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3.2 Perturbations 



To study the evolution of the perturbation, we combine the Eqs. ^ 
equation of motion of the gravitational potential 

$" + 2in - ^)$' + 2{n' - n^)^ - v^^ = 8ttG{2h 



TU|) and pT|) and obtain the 
^)V'<5V' . (20) 



Since there is only a single physical scalar field degree of freedom for adiabatic scalar metric fluctua- 
tions, all other perturbation variables can be determined from knowing the evolution of the background. 
A frequently used variable is the curvature fluctuation in comoving coordinates, which is given by 



n 



(21) 



In the case of an expanding universe with regular matter, this variable is known to well describe the 
adiabatic fluctuations on large scales since it is conserved on super-Hubble scales because it obeys the 
equation 

2fc2($' +H$) 



{l + w)C = 



However, this equation becomes singular both when the universe is at the bounce point when the Hubble 
parameter H approaches and also when the universe is crossing the phantom boundary w = — 1. 
Therefore, ( is ill-defined near the bounce point, as already remarked in [301 SI]- Thus, we will focus on 
the evolution of $ which is well defined throughout the bounce. We are focusing on adiabatic fluctuations 
only in this paper. Since there are two matter fields, it is possible to have entropy fiuctuations. We hope 
to return to this question in future work. 




Large k Small k 



Figure 3: A sketch of the evolution of perturbations with different comoving wave numbers k in the 
Quintom Bounce. 



Since in the linear theory of cosmological perturbations all Fourier modes evolve independently, we 
will follow each mode, labelled by its comoving wave number fc, independently. Interestingly, in our 
model there are four times when a perturbation mode crosses the Hubble radius. The evolution of scales 
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is sketched in Fig. [31 Perturbations start out inside the Hubble radius in the far past. Since in the initial 
heating phase the Hubble radius is shrinking faster than the scale factor, modes exit the Hubble radius. 
Once Period 2, the "slow climb phase" starts, the scale factor shrinks nearly exponentially whereas the 
Hubble radius is constant. Thus, wavelengths which are not too large will re-enter the Hubble radius 
during this phase. During the rest of this phase, and throughout the bounce period, the scale remains 
sub-Hubble, and exits the Hubble radius again in the post-bounce inflationary phase, only to re-enter the 
Hubble radius at late times. Modes with an extremely long wavelength will remain outside the Hubble 
radius throughout Period 2. Since the Hubble radius diverges at the bounce point, even these scales will 
enter the Hubble radius before the bounce. 

To summarize, scales k fall into two broad categories, as sketched in Fig. [31 The flrst consists of 
modes with comoving wave number k large enough for the perturbations to return inside the Hubble 
radius before the bouncing phase and escape outside after the bouncing phase; the second consists of 
modes with small k such that they enter the Hubble radius during the bouncing phase. We will show 
later that there is a difference in the post-bounce spectra of fluctuations between these two range of scales 
these. 

In a bouncing universe, it is unsuitable to use the equation of motion for Mukhanov-Sasaki variable 
V [m [33] (related closely to () because this variable is singular both at the bounce point and the 
cosmological constant boundarv [iOl 1^ 1 1 8j . Therefore, in the following we will directly study the behavior 
of perturbations from Eq. ([^ . 

In the following we will develop analytical approximations to study the evolution of the fluctuations 
in the various phases. We will denote the time corresponding to the end of inflation by tc, the end of the 
bounce phase by ts- 

3.2.1 Heating Phase 

Long before the bounce, when t <C — ts, we have an average value of the EoS parameter which is w = 0, 
and the contribution of ip can be neglected. In this region, we can treat the evolution of the universe as 
being dominated by a single field (p, with the Friedmann equations being 

n'^^ilcP'^ + a'V), (22) 

n'o.^-^{-cb" + a'V), (23) 
using conformal time and Hubble parameter. Further we have 

r 2UU' - n" ^24) 



(f)' 2{rp - H') ' 

and a (X. rf' , with 77 < during the contracting phase. Thus, we obtain the expressions 

2 2 4 

Making use of these relations, we obtain the following equation for the gravitational potential 

$1! + + kHk = , (26) 

in momentum space. Note that this equation is very similar to the perturbation equation in a matter 
dominant period except that the last term k^^k does not exist for hydrodynamical perturbations (ordinary 
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hydrodynamical matter is non-relativistic, our scalar field matter is relativistic) . The analytical solution 
of the above equation of motion can be written as follows, 

= V'Hk-^ AkJiikj]) + k^^ BkJ^likT])] , (27) 

where the coefficients Ak and Bk can be determined from the initial condition. 
To set the initial conditions, we first define the variable 

This is the variable in terms of which the initial conditions are determined by canonical quantization |38| . 
We will take the Bunch-Davis vacuum as the initial conditions (when rj is sufficiently far away from the 
bounce point) 

Uk = -^e-^'^ . (29) 
V 2fc'' 

During Period 1 when w — we have ^ ~ 77^'^. Therefore the initial condition for the gravitational 
potential becomes 

<^,^Atj-'^ , A = 47rGV^773, (30) 

where A is a normalization factor. When |fc?7| ^ 1 we can take the asymptotical form of the Bessel 
function and, matching with the Bunch-Davies vacuum, determine the coefficients. The result is 

Ak - i^Aki , Bk = -^Ak-i . (31) 

When the wavelength of perturbation is larger than Hubble radius with |fc| <C 7i, then by using another 
asymptotic form of the Bessel functions we obtain the following form of ^k'. 

From this result, we can see that of the two modes of $ one is constant and the other growing during the 
contracting heating phase. 

3.2.2 Slovir-Climb-Contracting (SCC) Phase 

For the universe in the contracting phase, we have H < 0. In this case, 3i?0 is anti-frictional, and instead 
of damping the motion of (j) in the expanding phase it accelerates the motion of (j). Consequently, the 
canonical scalar field (j) will eventually stop oscillating and start to climb up its potential. We denote 
this transition time as r]c. For {rj > rj^) on, the universe is in a period of SCC phase with w ~ — 1. We 
assume that in the beginning of this period the ghost field can still be neglected. Therefore, the equation 
of motion for the gravitational potential becomes 

+ 2aH<i>'k + 2{a - e)n^^k + k^^k = , (33) 

where we take the slow-climb (slow- roll) parameters as 

H H 

e = - — r and a = . (34) 

2HH 



ttWe will discuss other possible initial conditions in a foUowup paper. 



11 



As is well known, we have 



(35) 



where a subscript 'c' indicates that the corresponding quantity is evaluated at the time rjc- 

Making use of the above relations, we can solve the equation p3p for the gravitational perturbation: 

<^u^{ri~ficT[k~''CkJ,{k{'q~Vc)) + k''DkJ-Ak{v~flc))] , (36) 
i and f]c = rjc + l/H-c- We can obtain the asymptotic form of Eq. ([36| in the 



where a ~ v 
sub-Hubble region 



2 C}^ 

-[— sin(fc(r/ - i)c)) + -Dfe cos(fc(r/ - i)c))] 

TT k 



and in the super-Hubble region 



(37) 



(38) 



In the latter form the two modes are a constant mode Dk and a mode with coefficient Ck which is growing. 

Through matching $ and C (these are the matching conditions derived in [3S]) in the heating 
phase and in the SCC phase at the time rjc, we can determine the coefficients Ck and as follows, 



Cfc 



-He 



^(l-^e)Afc + 3(l + e)i?,r;-5 



Dk = e{-Ak-3BkVcl 







(39) 



Therefore, in the SCC phase the main contribution to the perturbation comes from the initial growing 
mode of $fc. 



3.2.3 Bouncing Phase 

Once the universe enters the quasi-de-Sitter contracting phase, the contribution of the ghost field tp 
becomes more and more important, since its kinetic energy density scales as a^^. Therefore, the SCC 
phase will cease at some moment rjs- when the contribution of ip begins to dominate, and then the EoS 
parameter of the universe will cross w = — 1 and fall to negative infinity rapidly. Correspondingly, the 
energy density and the Hubble parameter reach zero when the bounce takes place at the time rjB ■ Since at 
this moment the Hubble radius 1/H approaches infinity, all wavelengths will be inside the Hubble radius. 
Thus, the perturbations will be once again oscillating close to the bounce point. Short wavelength modes 
will have time to perform many oscillations which are sub-Hubble near the bounce point, whereas long 
wavelength modes (which have a longer intrinsic period and also spend less time being sub-Hubble) will 
not have time to perform a complete oscillation. As we will see, this difference leads to a difference in 
the way the initial spectrum of fluctuations transfers through the bounce. 

It is rather complicated to solve the perturbation equation exactly from Eq. (|20p . In order to gain 
some analytical insight into the evolution of fluctuations in this region, we resort to some approximations 
in order to simplify the equation. 

A convenient parametrization of the Hubble parameter near the bounce point (we choose t ~ to 
correspond to the bounce point) is 

H{t) = at , (40) 

where a is a positive constant. 
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In this case, we obtain the analytic forms of the scale factor and the comoving Hubble parameter: 



where as denotes the scale factor at the bounce point. 

In the bouncing process, the field (j) is initially climbing up to the hill. However, as the absolute value 
of H is decreasing, (f) changes its sign which then leads to the field (j) reversing its direction and beginning 
to slowly roll down the potential. At this stage, the universe has bounced and H is positive. 

Making use of the parameterization H ^ at (see (HH])) we can solve Eq. Q to obtain 

i = -3at^y{^'-£^)^~3H , (42) 

around the bounce point. Consequently, we can derive another formula (j)" ~ — 27i(/)'. In Sec. 13.2.51 we 
will show that this is a very successful simulation by comparing with the numerical calculation. In order 
to show this point more explicitly, in Fig. H] we present a plot of the factor 3H + ^ in the bouncing 
phase. The values shown in the figure support the approximation we have made in Eq. (42) since from 
the figure, it can be seen that this factor vanishes around the bounce point. 




Figure 4: A plot of the evolution of the factor 3H + ^ in the bouncing phase. 



We suppose that the approximations above should be valid only in the neighborhood of the bounce 
point, and hence the square term and higher order terms of |?7 — 77_b| should be neglected. 
Substituting these approximations into Eq. (|20p . we obtain the following equation: 



where we define the parameter yi = ^aog- The solution to this equation can be written as 

= e-i^^^"-"-)' Uh,,,(VW2(^ - Vb)) + F, iFA-y, \; \y,{v ~ Vb?) 



(43) 



(44) 
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which is constructed by a m-th Hermitc polynomial and a confluent hyper-geometric function with m = 

1-2 
3yi 



- — ^ and two undetermined coefficients Ek, F^. These two functions are linearly independent, but both 



of them show oscillating behavior when yi/k"^ is sufficiently small. This is to be expected in our bounce 
phase since the Hubble radius blows up at the bounce point and all scales become sub-Hubble and hence 
the associated mode functions should oscillate, even for \k{'q — ris)] ^ 1- 

If the parameter ffi is big enough. We can make use of the following approximate 

$fc~e-32'i |i?fcC0s[fc(r/-77B)]+^fcSin[fc(r/-?7B)]} ^ (45) 

where 

o m+l rz 

Ek = , ^ - Ek (46) 

Fk^^^^E, + F,. (47) 

Interestingly, we obtain a factor e'l^^'''"''^^ which can enlarge the amplitude of the perturbation before 
the bounce and suppress it after the bounce. We sketch this solution in Fig. [5] and will find that it is 
consistent with the numerical results in Sec. 13.2.51 



and 




Figure 5: A sketch of the evolution of the gravitational potential $ in the bouncing phase. 

Note that there are two possible ways for the modes in the SCC phase to evolve into the bouncing 
phase, as described in the beginning of Sec. 13.21 (also see Fig. [3]). The first possibility is that the comoving 
wave number k is large enough such that the mode reenters the Hubble radius in the SCC phase. The 
second is that k is small enough and the mode re-enters the Hubble radius only during the bouncing 
phase. 

In the first case with k ^ Hb-, the spectrum has been oscillating already at late times in the SCC 
phase. Through matching <I> and C in the SCC phase and in the bouncing phase at the time rjs- we can 
determine the coefficients E^^^ and pj.^^ as follows, 

^(1) ^ ^iviivB- -vnf |(^ sm[k{r]B- - f]c)\ + E)k co^kij^s- - f/c)]) sin [fc (775- - rys)] 
^ ^ -[(Cfc+7iB_i?fe)cos[fc(77s--?7c)] 
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+ - kDk) MKVB- - f]c)]] cos[fc(r/B_ - rys)] I , (48) 

pm ^ ^iviim- -^nf |(^ sin[fc(ryB_ - 77,)] + co^kij^s- ~ ffc)]) cos[fc(r;s- - m)\ 

+Tj T^i^K^fe + ^B-^fe) cos[fc(77s- - flc)] 

k(r]B- - Vb) eriB- 

+ i^Ck - kDk) MKVB- - fic)]] MKVB- - 775)] I • (49) 

From this result we can read off that, in the case when k ^ TLb- the coefficients of the perturbations in 
bouncing phase depend in an UN-suppressed way on the coefficient Ck of the growing mode in the SCC 
phase. However, this is not an absolute relation since the coefficient Dk of the constant mode in the SCC 
phase may play an important role when we fine-tune the phases k{riB- ~ "Uc) and k{riB- — ^jb)- 

Secondly, we consider the case when k ^ "Hb-- Through matching the spectrum in the bouncing 
phase with that in the SCC phase, the coefficients can be determined as follows, 

^(2) ^ e3^'i(''--"-)'y||(i^fc-^)sin[fc(,7B--r7B)] 

■coa[k[rjB- - Vb)] 



k{r]B- - Vb 



-ehdvB-VB? VVvr ^1 cos[fc(77B_ - r^g)] , (50) 
k(riB- -ilBj e 



Pl^) ^ et^^(''--"-)'^|(i^fe-^)cos[fc(r;s--,7s)] 
+ 77 ^- -— sin[fc(?7B_ -11b)] 

^ ehi(vB-vnf^/N^Rlsm[kivB--VB)]. (51) 
k[VB- - VB) e 

From these two equations we find that it is usually the coefficient Dk of the sub-dominant mode in the 
SCC phase which mainly dominates the coefficients of the perturbations in the boimcing phase, unless 
we fine-tune the phase k[vB- — Vb)- 

3.2.4 Slow-Roll-Expanding Phase (SRE) 

Due to symmetry of the Quintom model, the scalar field (j) begins to slow roll down to the bottom of its 
potential and the ghost field -0 decays soon after the bounce. Therefore, a phase of slow-roll expansion 
(like slow-roll inflation) starts at the time vb+- We denote the scale factor at this moment as aB+- 
During this period, the Hubble parameter is nearly constant which is similar to what occurs in the 
usual inflationary scenario. Interestingly, the EoS of the universe evolves from below —1 to above which 
provides a link with the normal expanding history of our universe. 

The equation of motion for the gravitational potential in this phase is the same as in the SCC phase 
(Eq. ([55)) ). Thus, the solutions takes the form 

$fe = {V- VB+)"[k-''GkJ,Akiv - VB+)) + k'HkJ-AHn - VB+))] , (52) 
where vb+ = Vb+ + ^/'kiB+ and the values of a and 1/ are the same as those which appeared in Sec. 
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13.2.21 We can obtain the asymptotical form of Eq. ([5^ in the sub-Hubble region 

/2 G 
-[-^ sm{k{r] - fiB+)) + Hk cos{k{r] - fiB+))] (53) 
TT k 

and in the super-Hubble region 

$fe ~ -y/|[Gfc (77 - fiB+) + Hk] . (54) 

The second mode is constant (it corresponds to the usual dominant mode on super-Hubble scales in an 
expanding universe) and the first mode is decreasing. 

When determining the coefficients of the two perturbation modes, we again need to match with the 
perturbation in bouncing phase in two cases. In the case k ^ Ti.B+ and assuming that all the phases are 
nontrivial, we obtain 

^ e-f^i . - .^ k{EkSin[k{T]B+-VB)]+FkC0s[k{7]B+~VB)]} ,(55) 

2sin[k(r]B+ - riB+)\ J 

^(1) ^ ^-h^ivs^-vsf , sm[k{rjB+ - Vb)] + ffe cos[k{,jB+ - Vb)]} ■ (56) 

2cos[fc(r/B+ - r]B+)\ J 

Therefore, in this case both coefficients are of the same order which means that the two modes are both 
important in the sub-Hubble region. However, as soon as the perturbation evolves outside the Hubble 
radius, the decreasing mode becomes negligible. The amplitude of the dominant mode in the expanding 
phase is set by the amplitude of the growing mode in the contracting phase. We will compare our result 
with that obtained in other works on non-singular bouncing cosmologies in the Discussion section. 
In the other case we take k ^ Hb+ and also assume all the phases are nontrivial. We then obtain 

~ -e-iy'^'^^+-^^^"nB+^{EkSm[k{T^B+-VB)]+FkCos[ki7jB+-VB)]} , (57) 

Hi^^ ~ _ VB)e^k {Ek cos[fc(r7s+ - Vb)] - Fk sin[k{vB+ - rw)]} 

~ . (58) 

Thus, we obtain the result that for these very long wavelength mode, the perturbations contain, to leading 
order in k, only the decreasing mode at late time. This is similar to the behavior of fluctuations in models 
with a singular bounce, such as the Pre-Big-Bang [7] and Ekpyrotic universe |H] scenarios. In the Pre- 
Big-Bang scenario, it has been shown [Hj that the growing mode of the ffuctuations in the contracting 
phase matches (up to correction factors which are suppressed by powers of k) only to the decaying mode 
in the expanding phase. The amplitude of the dominant mode in the expanding phase is thus set by the 
amplitude of the constant (and thus sub-dominant) mode in the contracting phase. The same result also 
holds [m [321 1371 13H] in the case of four dimensional toy model descriptions of the Ekpyrotic scenario, 
although in the case of the Ekpyrotic scenario the physical model is defined in higher dimensions, and 
an analysis of the evolution of ffuctuations in the higher dimensional context [33 EU] shows that there is 
a non-trivial mixing between the growing mode in the contracting phase and the dominant mode in the 
expanding phase. 

3.2.5 Numerical Results 



In the above sections we have present the analytic calculation of our model, during the process we have 
taken an assumption that the right hand side of Eq. ([^0)) is negligible during the whole evolution. This is 
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a key to solving the perturbation equation in our model explicitly. As is pointed out in previous sections, 
the energy density of the phantom field is usually enough small in the heating phase, the SCC phase, the 
SRE phase and the cooling phase and hence we have -0' — 0. Therefore it is a good approximation to 
neglect the right hand side of Eq. ([20]) in these four phases. Moreover, we have argued that the relation 
2H+ ^ ~ in the bouncing phase which makes the right hand side of Eq. ([20]) negligible. To make sure 
the assumption valid during the whole evolution, we plot the ratio of the r.h.s to the term /c^$ appearing 
in the equation as shown in Fig. [6l From this figure we can find that the value of the ratio in our model 
is always much less than 1 which agrees with our assumption very well. 




Figure 6: Plot of the evolution of the absolute value of the ratio of the r.h.s to the term k^^ in Eq. ([20 
The initial values of the background parameters are the same as in Fig. [1] 



In order to make the previous analysis exact and precise, we also solved the perturbation equations 
numerically. We provide the numerical result in Fig. [3 Fig. [HI and Fig. [D In Fig. [71 we plot the 
absolute value of the Newtonian gravitational potential |<I>| for different values of k. We can see from this 
figure that in the period far before the bounce point, those curves are similar. At the beginning of the 
evolution, when the wavelength is sub-Hubble, the perturbation is near its Bunch-Davies vacuum and 
its value oscillates with increasing amplitude. One may notice that the initial oscillatory behaviors of 
different modes are the same as others which only depend on the background evolution in Fig. [7] We 
give the explanations as follows. In our model, the canonical field (j) initially oscillates around its vacuum 
in the contracting period which results in a heating phase. During this phase, the EoS oscillates between 
— 1 and 1 as shown in Fig. [21 In our analytical calculation we simply take the average value of the EoS 
10 = and show that there are two modes for the gravitational potential with one being constant and 
the other being growing mode. However, when we take the numerical calculation, we encounter the case 
that there has to be a point satisfying cf) — and w = —I in every oscillation. At this moment, there 
can indeed be a large growth of fiuctuations due to parametric resonance. Therefore it seems that there 
are oscillations on the evolution of the gravitational potential, but it is not true. Since these resonances 
are determined by the background evolution, those "oscillations" have the same frequency and features 
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regardless of the wave number k. By the same reason we can understand the periodical amphfications of 
Id as they appear in Fig. [9l This issue will be discussed in future work. After the mode exits the Hubble 
radius, the perturbation will be dominated by its growing mode, as is shown in Eq. ([5^ . When it enters 
into the SCC phase, it grows monotonically at first and becomes oscillatory near the bounce point after 
the wavelength re-enters the Hubble radius. 

During the bouncing phase, the curves of different fc-modes behave very differently. For large k, the 
oscillatory behavior starts earlier and lasts longer because the wavelength of such modes enters into the 
sub-Hubble region much earlier before the bounce and escapes from the Hubble radius much later after 
the bounce. The amplitude of the oscillations near the bounce point increases before the bounce and 
decreases after the bounce, which follows from the Hubble term in the equation of motion which acts 
as anti-friction in the contracting phase and friction in the expanding phase, as shown in Eq. (I45|) . On 
the other hand, the oscillating period for small fc-modes becomes much shorter since these modes do not 
enter in the sub-Hubble region until after the bouncing phase has started. 

Finally, after the bounce the universe enters into the SRE phase and the wavelengths of the fluctu- 
ations become super-Hubble once again. Now, there are also two very different evolution scenarios for 
fluctuations with different fc-modes. In the expanding phase, for large enough values of fc, the gravita- 
tional potential is dominated by the constant mode; while in the case of small values of fc, it is dominated 
by decreasing mode. 

We also plot the values of <i> of different modes around the bounce point in Fig. [SI We can see 
that, for larger values of fc, the oscillations begins much earlier with their amplitude being enlarged at 
the bounce point, whereas for smaller values of fc there is only a very short period of oscillation (with 
amplified amplitude) which is caused by bounce. 

After the bounce, the gravitational potential for large fc evolves almost as a constant. However, for 
small fc there is a decrease on the gravitational potential. All these features are consistent with the 
analytical results discussed above. 
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Figure 7: Plot of the evolution of the absolute value of the gravitational potential |$| in our model. The 
initial values of the background parameters are the same as in Fig. [TJ 
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Figure 8: Plot of the evolution of the gravitational potential $ in our model. The initial values of the 
background parameters are the same as in Fig. [T] 

Finally, we also provide a plot of the curvature perturbation ^ in Fig. From this figure we can see 
that there inevitably are oscillations of the perturbations near the bounce point. 



4 Discussion And Conclusions 

Bouncing non-singular cosmologies, with an initial contracting phase which lasts until to a non- vanishing 
minimal radius is reached and then smoothly transits into an expanding phase, provide a possible solution 
to the singularity problem of Standard Big Bang cosmology, a problem which is not cured by scalar-field- 
driven inflationary models. 

In this paper, we have considered the evolution of cosmological fluctuations in a particular bouncing 
cosmology, namely the Quintom Bounce. In this model, space-time continues to be described by Einstein's 
General Relativity. The violations of the null energy condition required to get a bounce are obtained 
by making use of Quintom matter, which allows a transition of the EoS w through the cosmological 
constant boundary. We have studied the evolution of the gravitational potential both analytically and 
numerically in a particular model of a Quintom Bounce in which matter consists of two scalar fields, one 
of them a Quintessence-like field, the other a Phantom-likc field. We have found that the evolution of the 
gravitational potential $ is completely non-singular. Whereas in models with a singular bounce such as 
the Pre-Big-Bang [71 and the Ekpyrotic fS^ scenarios, the perturbations must be matched across a singular 
space-like surface, which leads to serious conceptual problems [37], in our case no such problems arise 0. 
In an analytical approximation scheme, we can apply matching conditions at non-singular surfaces which 
are then well justified, and we can also study the evolution numerically without any approximations. 

Our second main result is that for short wavelength modes, modes which re-enter the Hubble radius 
in the slow contracting phase, the dominant fluctuation mode is determined by the growing mode in the 

W Recently, a ghost condensation mechanism to smooth out the singularity in the Ekpyrotic scenario has been introduced 
[511152) . and it has been shown in this context, using the techniques developed in |53II54|[55] that making use of iso-curvature 
modes, the dominant spectrum in the contracting phase can be passed on to the dominant spectrum in the expanding phase. 
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Figure 9: Plot of the evolution of the absolute value of the curvature perturbation |C| in our model. The 
initial values of the background parameters are the same as in Fig. [TJ 

contracting phase. This result is similar to what has recently been found in other non-singular bouncing 
cosmologies [301 SIl [53] . The first of these is a bounce in the context of mirage cosmology, the second in 
the context of a specific higher derivative gravity model, and the third in the context of K-essence. Earlier 
work on fluctuations in models with a kinetic term with opposite sign generating a nonsingular bounce is 
in [57] . The authors of [57] observed that the way the post-bounce spectrum depends on the pre-bounce 
spectrum depends sensitively on the details of the bounce. In this light, it is not surprising that other 
bounces show that the dominant mode in the pre-bounce phase does not couple to the dominant mode 
in the post-bounce phase [TU O [5111111 EOl EI] ■ For long wavelength modes, however, we find that the 
dominant fluctuation mode in the contracting phase at leading order in k only sources the decaying mode 
in the expanding phase, a result similar to what is obtained in singular bouncing models [^H51HSllT7lH5] . 

Thus, the fluctuations in the Quintom Bounce model evolve in way which combines aspects found in 
some other recently studied non-singular bouncing cosmologies [IHl HI] (the behavior of short wavelength 
modes) with aspects found in singular bouncing cosmologies (the behavior of the long wavelength modes). 
Which category of modes determines scales accessible to current cosmological observations depends on 
the detailed parameters of the model. Our model may help us to understand the physics of singular 
bounce and non-singular bounce more clearly. As a side remark, let us add that in the Quintom Bounce 
model there is no trans-Planckian problem [S21 ES] for fluctuations since, unless the periods of slow- 
climb contraction and slow-roll expansion last many Hubble expansion times, the wavelengths of all the 
observable perturbations will be larger than the Planck length at the bounce point and thus never enter 
the zone of trans-Planckian ignorance. 

Let us close with some more general comments related to the Quintom Bounce. Although this theory 
is still in the stage of development, there have been works done in the context of a bouncing scenario 
related to predictions for observations. For example, in [64] (see also [651 [66]) it was pointed out that with 
a bounce followed by slow-roll inflation it is possible to give a reasonable explanation for the suppression 
of the low multi-poles of the CMB anisotropics. In [T5] and fS7\ the possibility of obtaining a bounce in 
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a Quintom model with a single matter field with higher-derivative terms was explored. See also [5H1 El] 
for other works. 
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